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Abstract. Once the fractional quantum Hall (FQH) state for a finite size system is 
put on the surface of a cylinder, the distance between the two ends with open boundary 
conditions can be tuned as varying the aspect ratio 7 . It scales linearly as increasing 
the system size and therefore has a larger adjustable range than that on disk. The 
previous study of the quasi-hole tunneling amplitude on disk in Ref. [11 indicates that 
the tunneling amplitudes have a scaling behavior as a function of the tunneling distance 
and the scaling exponents are related to the scaling dimension and the charge of the 
transported quasiparticles. However, the scaling behaviors poorly due to the narrow 
range of the tunneling distance on disk. Here we systematically study the quasiparticle 
tunneling amplitudes of the Laughlin state in the cylinder geometry which shows a 
much better scaling behavior. Especially, there are some corssover behaviors at two 
length scales when the two open edges are close to each other. These lengths are also 
reflected in the bipartite entanglement and the electron Green’s function as either a 
singularity or a crossover. These two critical length scales of the edge-edge distance, 
Lz 1 and L£ 2 , are found to be related to the dimension reduction and back scattering 
point respectively. 
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1. Introduction 

The strongly correlated electron system reveals a plenty of non-trivial properties which 
beyond the single-particle picture. The fractional quantum Hall effect (FQHE) [1] 
is a paradigm of strongly correlated system that occurs in two-dimensional electron 
gas with a perpendicular magnetic field. The FQH state is one of the mostly studied 
object in the condensed matter physics which has topological protected ground state 
and non-trivial excitation. Especially, the FQH states on the second Landau level, 
such as v = 5/2 and v = 12/5, are expected to have non-Abelian excitations and 
have potential applications in topological quantum computation mm- Quasiparticle 
tunneling through narrow constrictions or point contacts that bring counter-propagating 
edges close could serve as a powerful tool for probing both the bulk topological order 
and the edge properties of fractional quantum Hall liquids. In particular, interference 
signatures from double point contact devices may reveal the statistical properties of 
the quasiparticles that tunnel through them [5], especially the non-Abelian ones P Cl¬ 
in disk geometry, a quasiparticle can tunnel from the center to the edge by a single 
particle tunneling potential Vtunnei = Vtd(d) which breaks the rotational symmetry j8j. 
The ring shape of the Landau basis wave function with angular momentum mh on 
disk, i.e., ~ r m e H r l 2 / 4 has radius \/2 mis. Therefore the tunneling distance 

d is tuned by inserting N qb flux quantum, or N q h quasiholes at the center, namely 
d/l B = \f2(N orb + N qh ) — ^2N qh where N orb is the number of orbitals. The shape of 
the system evolves from disk to annulus and finally to a quasi-ID ring as increasing 
N qbl . In the ring limit (or CFT limit) with d —)■ 0, or N qb —>■ oo, we found a universal 
analytical formula for the tunneling amplitudes of the bulk quasihole [9] and the edge 
excitations [10] . On the other hand, the quasihole tunneling amplitudes were found to 
have a scaling behavior as a function of the system size N e and the tunneling distance d. 
Interestingly, the fractional charge and the scaling dimension appears in the exponents 
of the scaling function [9]. However, if we look carefully at the data of the tunneling 
amplitude in disk geometry, the scaling function works not very well for small d which 
was treated as a finite size effect [8] due to a limited number of electrons can be handled 
in the numerical diagonalization. It is also hard to look into this region on disk since 
N q h —> oo while d —> 0. 

In this paper, we alternatively consider the physics properties of the FQH liquid in 
the cylinder geometry. The cylinder has advantages that the distance between the two 
ends is proportional to the system size and can be tuned from zero to infinity smoothly 
by varying the aspect ratio 7 = L y /L x , where L y is the circumference on the side with 
periodic boundary and L x is the length of finite cylinder with open boundaries. As 
a comparison, in the disk geometry, the tunneling distance has a maximum which is 
the radius of the system R = y/2N orb and d ~ N orb /^/2N qh for N qh N orb . This 
is very inconvenient when we want to look at the small d region since numerous of 
quasiholes, or flux will be inserted at the center of disk. On the other side, when 
7 —> 0, namely in the thin cylinder limit, two adjacent Landau orbitals have practically 
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zero overlap. In this case, the Hamiltonian is dominated by the electrostatic repulsion 
which contains the direct interaction ((p m (p n \V\ , Tn l Pm) and the exchange interaction 
(TmTn\V\TmTn)■ The ground state is generally called a charge density wave (CDW) 
state, or Tao-Thouless (TT) crystal state [lHOEj on torus with the electronic occupation 
pattern 1001001001 • • • in order to minimize the electrostatic repulsion energy. The wave 
function of the FQH state [13] can be obtained by diagonalizing the model Hamiltonian 
with hard-core interaction, or the Hamiltonian only with V\ ^ 0 in the language of the 
Haldane’s pseudopotential. The more interesting case is when 7 —>• 00 or L y —> 00 , to 
keep the total area of the surface L x L y = 2irl%N or b invariant, or keep the total penetrated 
flux invariant, L x then approaches to zero. It means the two counter propagating edges 
at the two ends of cylinder are coming close to each other and the system finally evolves 
into a one-dimensional system. In this case, because of the strongly overlap of all 
the Landau orbitals, the Gaussian factors of each Landau wavefunction are the same 
and can be erased by normalization. In this one-dimensional limit, the FQH wave 
function can be described by the Jack polynomials and therefore all the results are 
the same as that we did on disk in the ring limit. The Jack polynomial is one of the 
polynomial solutions for Calogero-Sutherland Hamiltonian HI which can describe the 
Read-Rezayi Z^-parafermion states with a negative parameter a and a root configuration 
(or partition). Jack polynomial is a powerful method in studying the FQHE as it can 
construct not only the model wave function for Read-Rezayi series P5l ESI HZ!, but 
also the low-lying excitations [T 8 l [19]. Another advantage for cylinder geometry is the 
computational convenient comparing with either the disk or sphere geometries which 
was discussed in the density matrix renormalization calculation [20]. In this paper, we 
reconsider the quasiparticle tunneling with cylinder geometry especially in the region 
of small tunneling distance. Here the quasiparticle can tunnel from the one edge to 
another as sketched in Fig. CD Thus the tunneling distance equals to the length L x of 
the system. We find a richer structure in this region and two characteristic length scales 
appears not only on the quasiparticle tunneling, but also in the wavefunction overlap, 
bipartite entanglement entropy and electron Green’s function. 

The rest of this paper is organized as follows. In section II, we consider the tunneling 
amplitude with varying the length of the finite cylinder for e/3 and 2e/3 quasiholes in 
Laughlin state. In section III, the bipartite entanglement entropy, both in orbital space 
and real space are discussed. The results of the electron Green’s functions are discussed 
in section IV and summaries and discussions in section V. 

2. Quasiparticle tunneling for Laughlin state 

For electrons on a cylinder with circumference L y in y direction in a magnetic field 
perpendicular to the surface, the single electron wave function in the lowest Landau 
level is: 

1 r ik v y r -\ ( x+ky ) 2 


^jif) = 1 3 ) 


( 1 ) 
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Figure 1. (Color online) The sketch of the quasiparticle tunneling from one side to 
another side of the cylinder. The tunneling path is given by a potential Hunnei = VtS(y) 
which has length L x . 

2ir 

in which k y = — j, j = 0 , ± 1 , ±2 • • • are the transitional momentum along y direction. 
Ly 

Here the magnetic length Zg = \Jhc/eB has been set to be unit. For a finite size 
system, the number of basis states or orbits, N orb , equals to the number of magnetic 
flux quantum penetrate from the surface. Each orbit occupies an area 2nl 2 B . Therefore, 
the length in x direction for a finite system is fixed with a given aspect ratio 7 , namely 

L x /Ib 

To study the quasiparticle tunneling of the Laughlin state at v = 1/3, a quasihole 
with charge e/3 or 2e/3 is put on one edge of the cylinder as shown in FigiU Here 
the model wavefunction for Laughlin state can be obtained by diagonalizing the model 
Hamiltonian with hard-core interaction, or just V\ 7 ^ 0 in the language of the Haldane’s 
pesudopotential. It can also be obtained by using the Jacks with the so called root 
configuration “1001001001 •••”. The quasihole state for e/3 and 2e/3 is just the 
translated states with one and two sites along x direction respectively. Or we can 
say a quasihole is inserted at the left edge of the finite cylinder which is represented as 
roots “01001001001 ••• ” and “001001001001 •• ■” for e/3 and 2e/3 respectively in the 
Jack polynomial description. A simple single particle tunneling potential 


Vtunnel = V t 8(y) 


is assumed. It describes a tunneling path along the x direction and therefore breaks the 
translational symmetry in y direction. Then the matrix element is (Zc|Ftunnei|ni), which 
is related to the tunneling of an electron from the single particle state | m) to state | k), 
is (set V t — 1) 


v p (k,m ) = (k\V tUDne i\m) = e 


-fr( m-kf 


= e y 


( 2 ) 


It is clearly that v p (k,m) = e L v m L v k 1/1 = g-GAs) 2 / 4 where d is the distance 
between the two Gaussians. The many-body tunneling operator can be written as 
the summation of ones for single particle r = £(?/,). Then we can calculate the 

tunneling amplitude for many-body wave function T = ('F f/ /j|r|'I / o) in which the |\k 0 ) and 
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Figure 2. (Color online) The tunneling amplitude T as a function of tunneling distance 
L x for e/3 (a) and 2e/3 (b) Laughlin quasihole for system with 4—10 electrons. The 
data from different systems collapse into a single curve with a scaling function in EqJ5] 
The exponents are aA 3 = 2/3 (c) and a e F = —l/3(d) for two types of quasiholes. 
The insert figure shows the enlarged rescaled data for small L x . Some crossover points 
are labeled by arrows. 


are the ground state and quasihole state wave function respectively. In this section, 
we just consider the tunneling amplitudes for e/3 and 2e/3 quasiholes in Laughlin 
state at v = 1/3. The matrix elements consist of contributions from the respective 
Slater-determinant components \mi ,..., m^r) G Tq and \ki ,..., k N ) G ^ q h- There are 
nonzero contributions only when the two sets rri \,..., m n and k±,... ,k]y are identical 
except for a single pair m' and k' with angular momentum difference k' — ml — N 
for e/3 and k' — ml = 2 N for 2e/3 where N is the number of electrons. Therefore, 

e/3 -T- N 2 2 e/3 —^-(2AT) 2 

we have v p = e L y and v p — e L y . The tunneling amplitude in the second 
quantization can be written as: 

T = (T^|r|T 0 ) = y>ifc 2 • ■ ■ k n \8{y i )C^C m \m 1 m 2 ■ ■■m n ). 

i 

From the Eq. [21 it is known that the tunneling amplitude decreases exponentially as 
increasing the tunneling distance which is proportional to \m — k\. The distance of 
the quasiparticle tunneling of mang-body state, or the length of the cylinder L x is 
determined by the size of the system and the aspect ratio 7 . For a IV-particle system 
at fixed filling factor, the L x can be tuned from 0 to 00 by changing the aspect ratio 7 . 
As that on the disk, a numerous of quasiholes were added at the center which makes 
the radius changes from y/2N orb to 0. And with a given L y , the distance between two 
single particle orbitals on cylinder is a constant which makes the tunneling distance 
scale proportional to the system size, i.e., d oc N on cylinder comparing with d oc N 1 ^ 2 
on disk. The linear relation guarantees a smooth change while varying 7 . 

Fig.E^a) and (b) show the tunneling amplitudes for e/3 and 2e/3 Laughlin quasihole 
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as a function of the tunneling distance L x . When L x —> oo, or 7 —> 0, the system is in 
a thin cylinder limit and the ground state is a crystal-like state in which electrons are 
separated, then the quasiparticle can not tunnel from one side to another, i.e., T —> 0 . 
On the other side, when L x —y 0, or 7 —> 00 , all the single particle orbitals collapse onto 
each other which corresponding to the ring limit on disk, or the CFT limit in which case 
the geometry factor of the many-body wave function can be neglected. Our previous 
studies Um show that the tunneling amplitude for e/3 and 2e/3 quasihole in the CFT 
limit for a system with N electrons can be exactly represented as: 
re/3 = 17(1001001...01001) 

1 J 17(01001001... 01001 ) 


Q 


r 2 e/ 3 ( i y) = 2\N- 


Q 


001001001...001 
1001001001...001 

01001001001...001 
001001001001...001 


(3) 


where 17(1001001...01001) = 1 x 3 x 6 x ... and 17 | ^ ] = 17(//) 17(A). For example, 


2x3 


= 0.6. The numbers in the 


r e/3 (2) = 2igf = 1.5 and 1^(2) = 2 x 2 - x2x4x!i 
fractions are the position of the Is in Eq. [3j Formally, the tunneling amplitude for 
e/3 quasihole for Laughlin state has an algebraic expression r e / 3 (lV) = jj ), 

in which M — \/v — 3 for Laughlin state and the beta function B is defined as 
B(x,j3) = T{x)T(P)/r{x + P). Fig. |2](a) and (b) show that the tunneling amplitudes 
saturates exactly at these CFT limit values when L x —y 0. In the medium region of L x , 
the tunneling amplitude has a dramatical change from these CFT values to zero. The 
state in this region is close to the Laughlin state, thus the signal of the decreasing of the 
quasihole tunneling amplitude can be seemed as a measurement of a phase transition 
(PT)-like from the thin cylinder state with zero tunneling amplitude to the CFT limit 
with a finite tunneling amplitude. Here we should note that we use the terminology 
PT-like instead of PT since there is actually no phase transition in the ground state 
while varying L x . The topological properties of the ground state in CFT limit are the 
same as that in the thin cylinder limit [21, [221122] • As shown in Fig. [2j c) and (d), 
the data for different system sizes collapse into each other after the following scaling 
conjecture is applied 

T q (N, L x ) = T 0 N~ aq e( qL * /2elB)2 . (4) 


The exponent a q is related to the scaling dimension of the quasiparticle as a q = l — 2A q . 
The u — 1/3 Laughlin quasihole operator can be written as if q h = e l< ^ m /v3 with a 
primary charge bosonic held p in CFT. Therefore, the scaling dimension for e/3 and 
2e/3 quasiholes are A 6 / 3 = 1/6 and A 2e / 3 = 2/3 respectively and then cc 6 / 3 = 2/3 
and ct 2e / 3 = —1/3. In disk geometry | 8 j, the best scaling parameter for 2e/3 was 
a 26 / 3 = —0.4 which has a large deviation from the theoretical prediction. We think this 
deviation should come from the insufficient tunneling distance which is maximized at 
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Figure 3. (Color online) (a)Tlie normalized tunneling amplitude as a function of 
tunneling distance for e/3 (a) and 2e/3 quasihole (b). The insert figures are their first 
order derivations which have peaks at L x 2 ~ 4.2/^ both for e/3 (a) and 2e/3(b) by 
extrapolating to N —> oo. 


the radius of the disk. On the other hand, as we discussed above, the tunnel amplitude 
data near d — 0 is missed due to the huge number of quasiholes needed to be insert at 
the center. On cylinder, as showed in Fig. [2(c) and (d), the scaling conjecture of Eq[4] 
works perfect when L x is large than a specific value. The reason we are saying this is 
that if we enlarge the rescaled data in the small L x region as shown in the inserted figure, 
crossover behaviors from different systems occur around Lff ~ 2.5 Ib and L C2 ~ 4.2 Ib as 
shown by arrows in Fig. [21(c) . The crossover at L c f ~ 2.5 Ib is also remained the same in 
Fig. [2(d) for 2e/3 quasihole. Here we should note that there is no crossover for larger 
L x in Fig. [2(d). However, we observe that the scaling behavior is starting to be broken 

c' 

down near L x 2 ~ 5.2 Ib- The first crossover at L c f ~ 2.5 Ib can be explained a transition 
from two dimensional system to one dimensional system. The ID system corresponds to 
Calogero-Sutherland model |2Tj which actually is the origin of the holomorphic part of 
the FQH wave function, or the Jack polynomials Gang nz]. The dimension reduction 
of the FQH state was also considered in the composite ferimion systems [25, 26] . Or 
we can say that the system is in the CFT limit while L x < Lf. The second critical 
value L C2 ~ 4.2 Ib for e/3 and L x 2 ~ 5.2 Ib for 2e/3 is the transition point where the 
scaling behavior is broken down. This can be explained by the broken down of tunneling 
behavior between two independent edges due to gluing the two anti-propagating edges 
together when varying 7 . Or we can say that the L C2 is the length scale at which 
two edges start to interact with each other. When L x < L x 2 , there are back scatterings 
between the two anti-propagating edges. The different values of L C2 between the e/3 and 
2e/3 quasiholes should be from the size difference of them. Another way to extrapolate 
the critial point L C2 is renormalizing the data by its CFT value from Eq. [3] as shown 
in Fig. [3] Interestingly, the data for different sizes have a scaling-like behavior with 
collapsing onto one curve. The insert plots in Fig. [3] are the first order deviations of the 
normalized tunneling amplitudes. Again the first deviations have peaks at L C2 ~ 4.2 Ib 
both for e/3 and 2e/3 by extrapolating to the thermodynamic limit with N —» 00 . 
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Figure 4. (Color online)The ratio T 2e / 3 /T 2e / 3 as a function of L x . The solid line is 
plotted using Eqj5]and he dash line is the fitted line with T 2e / 3 /r e / 3 « e _0 ' 078 ( x / iB ) . 
The insert figure is the enlarged part with smaller L x for N = 10 electrons. The 
numerical data deviates from the asymptotic curve near L x 2 ~ 4.3 Ib- 


On the other hand, besides the N dependence of the tunneling amplitudes, the 
Eqj4] tells us the ratio of the two types of tunneling amplitude is expected to has an 
asymptotic behavior which depends on L x : 


p2e/3 

r e / 3 


_ ^—[(2L x /3) 2 —{L x /3) 2 ]/(2l B ) 2 


e -0.083 (L x /l B ) 2 


(5) 


In FigJU we plot the ratio r 2e / 3 /T e//3 as a function of L x for system with 6 — 10 electrons. 
Unlike the data on disk | 8 ] in which there was a sudden change while the first quasihole 
was inserted at the center, the ratio on the cylinder is smooth as a function of L x . 
The data in FigQ] can be fitted by a solid line with r 2 e// 3 /T e / 3 ps e -o.o78 (x/i B ) 2 w hich 
is consistent with the expected behavior in Eqj5] With comparison to the disk, the 
ratio of the tunneling amplitude for e/4 and e/2 for Moore-Read state on disk has 
an asymptotic T e / 2 /r e / 4 ~ e -o.o83 {d/i B ) 2 w hich has a relative large deviation from the 
expected behavior r e / 2 /T e//4 ps e^ 0 , 047 (d/i B ) 2 . And for Laughlin state, the numerical 
and theoretical prediction are r 2 e / 3 /r e / 3 ~ e -om(d/i B ) 2 ^ p 2 e/ 3 ^pe /3 ^ e -o.o83 (d/i B ) 2 
respectively. Moreover, if we just plot the data for 10 electrons as insert plot in Figjdl 
it is shown that the deviation of the asymptotic behavior occurs near L x = 4.3 Ib which 
is close to ~ L C2 . This deviation also demonstrates that the tunneling has been affected 
by the edge-edge interaction at this length scale. 


3. Bipartite Entanglement entropy 

The idea that the quantum entanglement [2?] in a bipartite system can describe different 
phases of matter has emerged over the past years. This approach has provided plenty 
of new insights, while traditional methods based on symmetry breaking and local order 
parameters in Landau theory are fail. More precisely, a bipartition of the quantum 
system is defined when the Hilbert space factors into two parts TL = FLa ®FLb- The 
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Figure 5. The entanglement entropy as a function of L x in OC (a) and RC (b) and 
as a function of L y in OC (c) and RC (d). The insert plot in (a) and (b) are the 
scaling of the position of the peak in the entropy. The (c) has even-odd effect in large 
L y region and the two fitting straight lines are 0.139 L y — 0.498 ± 0.001 (even) and 
0.155L y — 0.580 ± 0.004 (odd) respectively. The (d) has even-odd effect in small L y 
region and the fitting line for large L y is 0.198L y — 0.546 ± 0.006. 


bipartite FQH system can be implemented in both the momentum space and the real 
space of the two dimensional electron system. The former is called the orbital cut 
(OC) [28] and the later real space cut (RC) [29]. With a bipartition, a pure quantum 
state \ip) can be expressed in the form of Schmidt decomposition 

W) = Y e _6/2 | iff) ® |^f >, (6) 

i 

where | iff) and | iff) are orthonormal sets in Ha and Hb respectively and the value 
of fi in Schmidt singular values e - ^/ 2 are the entanglement “energies” in entanglement 
spectrum [30j. Equivalently, the reduced density matrix pa = tr j B|-0)('0| has eigenvalues 
A i = e~^\ The Von Neumann entropy 

S A = -tr[p A lnp A ] = - Y -MnAj = Y& e ~^ ( 7 ) 

i i 

generally scales linearly with the area of the cut between parts A and B and with a 
universal order 0(1) correction, namely the topological entanglement entropy [31., [32], 
[33] . i.e., S = aL — The topological entanglement entropy 7 * of the ground state 
for a fully gapped Hamiltonian is one robust measure of quantum entanglement in a 
topological phase in two dimension system. In the FQH state, 7 * = logV where T> > 1 is 
the total quantum dimension of the system. For Laughlin state at v = 1/3, the quantum 
dimension is V = \/3 and therefore 74 ~ 0.549306. 

On the cylinder, we intend to divide the system into two equal subsystem with the 
same number of orbitals in OC or the same length in RC. However, the number of orbitals 
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for IV-electron Laughlin state is N or b = 3N — 2 which has the same parity as N. Then 
the bipartition of the orbitals should has an even-odd effect, which can be defined as the 
orbital difference between Ha and Hb, namely \N° rb —Nff rb \ is 0 for even N and, 1 for odd 
N. Intuitively, the effect of orbital difference should be diminished as increasing L x due 
to the local properties of the entanglement entropy, or inversely, it becomes more clear 
in small L x , or large L y region. This can be shown in FigJS^c). In the RC case, it is easy 
to comprehend that even-odd effect exists, especially in the thin cylinder limit. Taking 
the N = 2 and N — 3 thin cylinder crystal-like states as examples, their wavefunction 
are single Slater-determinant ty TT (2) = 11001) and T tt ( 3) = 11001001) respectively. 
Then the position of the RC cut for state with even N has zero electron density which 
induce a zero entanglement entropy. On the other hand, there is an electron locates at 
the position of the RC cut for odd electron TT state. Then the electron density reaches 
its maximum and the entanglement entropy saturates at a specific value, which is the 
same as that for cutting a single electron wave function into two equal parts which is the 
classical Von Neumann entropy Stt = log(2) ~ 0.693147 which is shown in Figj5])d). 
When the cylinder is bipartite along the y direction, the L y is actually the length of 
the cutting, or the “area” between two subsystems. Then the topological entanglement 
entropy can be extrapolated from FigJHfc) and (d). We found in the case of the OC, 
the data for systems have the same parity are sitting on the same curve as a function of 
L y . For large L y , they can linearly be fitted and the topological entanglement entropy 
for even and odd parities are 7 * ~ 0.498 and ~ 0.580 respectively. The exact value 
7 1 = log \/3 is in between them. A more accurate extrapolation can be obtained in 
RC entanglement entropy which shown in Figj5])d) where the ~ 0.546 which is close 
to the exact value. An interesting phenomenon is both the OC and RC entanglement 
entropies saturate at zero (log(2) for odd parity in RC) near L y ~ L x 2 . Therefore, we 
conclude that L C2 is the length scale in y direction that the CDW behavior appears. 

In Figj5])a) and (b), we plot the bipartite entanglement with OC and RC as 
a function of L x . The OC entanglement entropy has a peak at L x ~ 4.6 Ib in 
the thermodynamic limit as shown in the insert plot. In the RC, the peak of the 
entanglement entropy is not as sharp as that in the OC since the entropy while 
L x — y 0 decreases very slowly. The position of the peak in the thermodynamic limit is 
L x ~ 3.7 ± 0.1 Ib- The error bar origins from the strong even-odd effect in this region. 
The difference between OC and RC can be explained by the different width of the cuts 
in realspace. The OC has a wider cut range and is more sensitive to the change of 
L x . It is known that the entanglement entropy has a singularity at the critical point 
of the QPT due to the divergence of the quantum fluctuation. However, since there 
is no phase transition occurs while varying the L x |2TJ 122) 123], the increment of the 
entanglement entropy origins from the correlations between two edges. Therefore, the 
two length scales L x ~ 4.6 Ib and L x ~ 3.7 Ib for OC and RC respectively should be 
related to the L° 2 ~ 4.2 l B . 
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Figure 6. (Color online)The log-log plot of the electron Green’s function \G(L X )\ for 
the Laughlin state for 5-12 electrons as a function of L x . The data in the large and 
small L x region can be fitted by ~ Lx 20 ' 8 and ~ L x respectively. L x = 2.5 Ib is the 
saddle point between these two phases. 


4. Electron Green’s Function 


Tunneling characteristic at the edge has long been regarded as an experimental method 
of measuring the topological order of the FQH liquids. For tunneling from a three- 
dimensional Fermi liquid into the FQH edge, chiral Luttinger liquid theory [33j leads 
to a non-Ohmic tunneling I — V relation / oc V a with a ^ 1, in sharp contrast to the 
Ohmic prediction of a Fermi-liquid-dominated edge with a — 1. The electron Green’s 
function is defined as 


G( r - r') 


(l/’l^(r)d' e (r / )|?Q 

(?/#) 


( 8 ) 


where the T+(r) and T e (r / ) are held operators which create and annihilate an electron 
at position r and r' respectively. If we consider the tunneling path |r — r'| along the edge 
of the FQH droplet, the edge Green’s function shows a scaling behavior with a — 1/v 
for long distance tunneling [ 351 [36[| . 

In this section, we consider the electron tunneling from the left edge of the cylinder 
to the right one, namely the electron correlation function between two anti-propagating 
edges as a function of L x . The results are shown in Figjb] It shows that the Green’s 
function decreases dramatically when L x is larger than the saddle point which is the one 
dimensional limit threshold value Lf = 2.5 Ib- The data for large tunneling distance 
near L x ~ 10/^ obeys a power law behavior with an exponent less than —20. In 
the large L x limit, obviously, the Green’s function is zero in the Tao-Thouless state 
which is an insulator. We also checked the v — 1/5 Laughlin state and found that 
the electron Green’s function has the same power law behavior in this region. Thus 
we think that the exponent in the large L x region depends on the interaction between 
electrons. The electron Green’s function scales as a G(L X ) oc L x which has a positive 
exponent one. Generally, the electron Green’s function at zero temperature decays as 
G(r) rvj 7* 1 - Q [371138] in which a > 0 and a/2 is the “anomalous dimension of the 
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fermion operators. The case for a = 0 corresponds to the normal Fermi liquid and 
a > 0 is due to the correction of the electron-electron interaction which is the behavior 
of a Luttinger liquid. On the other hand, when L x < L c f , the system enters into a one 
dimensional phase which is described by the Calogero-Sutherland model. The reason 
that the correlation decreases as reducing L x is due to the repulsion between electrons, 
or strictly speaking, the electron Green’s function drop to zero in the ID limit. 

5. Summary and discussion 

As a conclusion, we confirm that the quasihole tunneling amplitude in the cylinder 
geometry obeys the scaling conjecture in Eqj4] and the scaling behavior is much better 
than that on disk. Generally the scaling behavior works well when L x > L C2 where 
L C2 ~ 4.2 Ib for e/3 and L C2 ~ 5.2 Ib for 2e/3 with a difference due to different size 
of the quasihole. The L C2 can be explained as the threshold value of the edge-edge 
back scattering between two edges. It appears not only in the quasihole tunneling 
amplitude calculations, but also in bipartite entanglement entropy. Therefore, the L C2 
is the smallest length scale that guarantees there are two independent edges at two 
ends of the cylinder. It should be the benchmark of the sample size in designing of the 
experimental setup of the quasiparticle tunneling and interference [39^ HQ]- Moreover, 
we found another critical value Lf ~ 2.5 Ib, which is universal for different types of 
quasiholes. It can be explained as the critical width evolving from a 2D system to ID 
system which is described by the Calogero-Sutherland model. Bipartite entanglement 
entropy has a singular behavior near L C2 due to a contribution of the edge-edge back 
scatterings. The topological entanglement entropy is extracted from the OC and RC 
entanglement entropies as a function of L y in a finite size system. The Lff plays a role 
of a saddle point in the single particle Green’s function where the system enters into a 
one dimensional description. The scaling exponent of the Green’s function is one while 
approaching to the ID limit. We notice that the Lf is actually the correlation length 
of the Laughlin state as mentioned in the iDMRG calculation [4l|. Here we should 
admit that we just consider the Laughlin state of the model Hamiltonian with hard¬ 
core interaction, or with V\ pesudopotential. For a realistic coulomb interaction or the 
FQH state in higher Landau levels, we believe that the similar behaviors exists which 
may at most have small modifications on the value of these lengh scales. 

We thank X. Wan and G. T. Liu for helpful discussions. This work was supported 
by NSFC Project No. 1127403 and Fundamental Research Funds for the Central 
Universities No. CQDXWL-2014-Z006. NJ is also supported by Chongqing Graduate 
Student Research Innovation Project No. CYB14033. 
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